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SUMMARY 

The  second  method  used  by  Lyapunov  to  stuiy  the  staolllty 
of  the  trivial  solution  of  dx/dt  ■  Ax  -f  f(x)  leads  to  tne  problem 
of  solving  the  matrix  equation  AX  XA'  »  C.  It  Is  shown  that 
this  question  Is  related  to  Kronecker  products  and  Kronecker 
sums,  and  the  general  problem  of  solving  AX  XB  »  C  Is  resolved 


In  this  fashion. 
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KRONECKHR  PRODUCTS  AND  THE  SECOND  METHOD  OF  LYAPUNOV 

By 

Richard  Bell^ian 

§1 .  Introduction 

In  the  second  metnod  ured  by  L:>ap.uno7  In  the  Inver  t  legation 
of  the  stability  of  tne  trivial  solution  of  tne  nonlinear  vector- 
matrix  system 

(1)  3*  ,  Ax  f-  f(x), 
dt 

the  question  arises  of  solving  the  matrix  equation 

(2)  AX  +  XA'  =  C, 

for  the  unknown  symmetric  matrix  X.  This  problem  was  treated  by 
W.  Hahn  In  a  recent  publication,  [_?]  ,  using  the  reduction  of  a 
matrix  to  triangular  form. 

The  same  problem  was  encountered  by  the  author,  [IJ,  in 

connection  with  the  evaluation  of  the  Integral  J  ■  ^/  )y,Dy)dt, 

0 

where  y  Is  the  solution  of  dy/dt  ■  Ay,  y(0)  ■  c,  and  A  is  a 
stability  matrix,  l.e.  all  cnaracterlstlc  roots  have  negative 
real  part.  The  question  of  existence  and  uniqueness  of  solution 
was  resolved  In  a  non— algebraic  fasr.lon  by  noting  that  the 
solution  of  (2)  Is  given  by 

CD  At  A  '  t 

(?)  X  -  -  e  Ce  dl, 

o 

under  the  assumption  tliat  A  Is  u  stability  matrix. 
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In  this  paper  we  shall  present  another  approach  to  this 
problem  using  the  concept  of  the  Kronecker  product  of  matrices. 
As  we  shall  see,  this  methoa  enables  us  to  resolve  the  general 
problem  of  determining  when  the  equation 

(4)  AX  +  XB  -  C 

has  a  solution. 

§2.  The  Kronecker  Product 

Let  A  and  B  be  two  nxn  matrices,  x  and  y  characteristic 
vectors  of  A  and  B  respectively,  and  A  and  ^  characteristic 
roots , 

(1)  Ax  ■  Ax,  By  ■  uy . 

2 

Let  z  be  the  n  -dimensional  vector  formed  as  follows 


r-1097 

-> 


Referring  to  (1),  It  Is  easy  to  see  that  z  lo  a  charactei Istlc 

2 

vector  of  tne  n  -dlmenflonal  matrix 


(3) 


A  X  B  » 


This  we  recognize  as  tiie  Kronecxer  product  of  tne  matrices  A 

•  2  ^ 
and  B.  Its  n  characteristic  roots  are  the  comblnatlona  , » 

with  tne  characteristic  vectors  obtained  as  above. 


§3*  The  Kronecker  Sum 

Let  Aj^  ■  I  -»-  €  A,  B^  »  I  €  B.  Then  the  Kronecker 
product  is 


(1) 


(I  -»■  €  A)  X  (I  ^  €  D) 


I  X  I  t  fell  X  B  4-  A  X  ij  .  0(6“^) 


(l  ■♦■6b^^)(l  4-€A)  ... 

(I  -^€b^^)(l  ^€A)  ... 


I  ^  €  (b^^  I  A)  .  .  .  \ 

I  €{b^^  I  >  A) 

:  / 


The  n— dimensional  matrix 


•  Kronecker  products  arise  also  in  the  study  of  the  product 
of  stochastic  matrices.  See  H.  Eellman,  Limit  Theorems  for  N'on- 
Commutatlve  Operations,  Duke  Math.  Jour.  Vol.  (1934),  pj  .  4;jl 


300 . 
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(2) 


A  0B  - 


A  +  b  ^  I 

A  -»■ 


A  +  bj^^I  .  .  A  +  b^^I 


we  ma>  call  the  Kronecker  sum.  Its  characteristic  roots  are 
and  Its  characteristic  vector?  are  as  given  In  (2.2). 


§4  .  ^  Alternate  Derivation 

This  matrix  may  also  be  obtained  by  considering  the  two 
linear  systems 

dx-  N  dy.  N 

(l)  ■  L  ax,,  «  Z  ^  “  1>2,...,N, 

dt  j-1  dt  j-l  > 


2 

and  forming  the  n  —dimensional  linear  system  satisfied  by  the 
vector  In  (2.2).  Since 


(2) 


d  ^^1 

— 

1  J  dt  J 


dt 


1 


dt 


we  obtain  the  matrix  In  (3*2). 


§5.  AX  f  XB  -  C 

Let  A  and  B  be  given  n— d Imens lonal  matrices,  and  consider 
the  matrix  equation 


(1)  AX  XB  -  C, 

where  X  Is  the  unknown  matrix. 

This  system  of  equations  for  the  unknown  elements  x,  ,  ol  X 

1  j 


has  the  form 


F-109'7 

^29-37 


It  follows  that  the  equation  In  (1)  possesses  a  unique 
solution  for  all  C  If  and  only  If 

(3)  ^^+yu,40,  l,J-l,2,...,n, 

where  the  are  the  characteristic  roots  of  A  and  the  the 
characteristic  roots  of  B. 

§6.  AX  XA'  =  C 

It  Is  easy  to  see  that  the  same  techniques  apply  to  the 
equation  In  (1.2),  so  that  a  necessary  and  sufficient  condition 
that  this  equation  has  a  unique  solution  for  all  C  Is  that 
^1  4-  Aj  t  0,  1,J  »  1,2, ...,n.  In  particular,  this  condition  Is 
satisfied  If  A  Is  a  stability  matrix. 
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